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Abstract. The group G is called n- rewritable for n > 1, if for each sequence of 
n elements x\,X2, ■ ■ ■ ,Xn G G there exists a non-identity permutation cr G Sn 
such that xiX2 ■ ■ ■ Xn = x„(i\x^r2) ' ' ' ^tT(n)- Using computers, Blyth and 
Robinson (1990) verified that the alternating group A^ is 8-rewritable. We 
report on an independent verification of this statement using the computa- 
tional algebra system GAP, and compare the performance of our sequential 
and parallel code with the original one. 
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rS^ • Let n > 1 be an integer. Following |1,, a group G is said to be totally n-rewritable, 

j^ I or have the rewriting property P„, if for each sequence of n elements xi,X2, ■ ■ ■ ,Xn 

of the group G there exists a non-identity permutation a ^ Sn such that 

xiX2---Xn = a;<T(i)a;o.(2) • • ■ x^(^„) ■ 

^ . Clearly, all abelian groups satisfy P2, and if G satisfies Pk then it also satisfies Pk+i- 

l/^ ' On the problem session of the conference "Arithmetic of Group Rings and Re- 

lated Objects" (Aachen, Germany, March 22-26, 2010) Eli Aljadeff (Technion, 
r^ ' Haifa, Israel) suggested the following problem: 

10 ' Prove that the alternating group A^ has the property Ps 

o ; 

^— *; ' He referred to the computer verification of this statement reported in [T], and 

demonstrated how to show that Ar, has the property Pio using group rings tech- 
nique. He also suggested that since [I] appeared twenty years ago, nowadays this 
result probably could be verified much faster. Motivated by this, the author veri- 
fied that ^5 has the property Pg using the computational algebra system GAP [2] 
5^ I and compared the performance of the sequential and parallel GAP implementations 

with the one described in [1]. 

To check that the group G is n-rewritable using the brute force approach, one 
may enumerate all n-tuples of distinct elements of G and check that each of them 
may be rewritten. Of course, even for A5 the number of tuples to check will be 
enormous, so this approach will not work. 

There is, however, a simple observation that allows to reduce the number of 
checks substantially. The algorithm described by Blyth and Robinson in [T] con- 
structs all non-rewritable words of length 2 which are non-equivalent with respect 
to the action of Aut(G'). On the next step, these words are used to construct all 
non-equivalent non-rewritable words of length 3, and then the process is repeated 
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until the first n for which there are no non-rewritable words of length n will be 
found. 

The next table contains information about the number N{r) of non- rewritable 
words of length r in A5, determined using computers and listed in [I]: 



r 


N{r) 


2 


29 


3 


1315 


4 


43121 


5 


528069 


6 


187719 


7 


1320 


8 






The authors of [1 wrote that first these data were computed over a period of two 
weeks by a PASCAL program on a Micro VAX II, and then verified by a parallel 
C-\ — h implementation that produced the same result on four Sun 3/60 machines in 
less than three hours. 

We were interested to compare the reported performance with the runtime that 
can be achieved nowadays in a sequential version of GAP on modern computers. 
To test our implementation, we used an 8-core Intel server, with dual quad-core 
Intel Xeon 5570 2.93GHz / RAM 48 GB / CentOS Linux 5.3. First we present the 
algorithm in the pseudocode, following its textual description from |T]: 



RewritabilityLength(G, m) 
1 A-h- AutomorphismGroup[G] 



2 


X -s- No 


ntrivialOrbitRepresentatives [A, G] 


3 


n^ 1 




4 


repeat 




5 




n <— n + 1 


6 




nrw <- EMPTY LIST 


7 




for u in X [> u is a non-rewritable word of length n — 1 


8 




do K ^ Intersection[ 

Stabiliser[A, u[l]], 

Stabiliser[74, u[n - 1]]] 


9 




if Size[X] > 1 


10 




then y ^ NONTRIVIALORBITREPRESENTATIVES[iir, G] 


11 




else y^G\{lG} 


12 




for V in y 


13 




do t = Concatenation [u,z)] 


14 




if ISREWRITABLEWORD[t] 


15 




then Append [nrw, i] 


16 




if nrw = EMPTY LIST 


17 




then return n 


18 




else X •<— nrw 


19 


until n = m 


20 


return 


fail 
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The pseudocode above refers to the fohowing procedures: 

• AutomorphismGroup[G] returns the automorphism group of G; 

• NontrivialOrbitRepresentatives[A, G] return the hst of representa- 
tives of orbits of non-identity elements of the group G under the action of 
its automorphism group A] 

• STABiLiSER[yl, g] returns the stabihser of an element g e G in the group A\ 

• IsREWRiTABLEWORD[t] checks if the word t is rewritable. 
Other names of procedures should be self-explanatory. 

The first implementation looked very much like the pseudocode above, and it 
took almost 34 hours to run (though it used one CPU, other CPUs were used 
for other jobs, so we can not guarantee exact measurement). The second version 
was optimised to achieve more efficiency on the stage when most of non-rewritable 
words of length k can not be extended to non-rewritable words of length k -\- 1. 
Concatenation of lists was replaced by changing the last element "in place" , and 
IsRewritableWord was insered directly into the loop without a call to a separate 
function that, in it turn, used ForAny. Additionally, intersection of stabilisers was 
computed in a loop which breaks if a trivial subgroup is constructed. This and 
some other minor optimisations allowed to reduce the runtime to about 15 hours. 
Finally, we traded space vs time and stored not only tuples, but also stabilisers of 
their elements in Aut(G). This permitted further speedup and reduced the runtime 
to be less than ten hours. At this stage we performed six clean measurements on 
a machine not running other user's jobs, and the average runtime was 9 hours and 
41 minute. 

The GAP code for the function RewritabilityLengthis given in the Appendix. 
As you can see from the example of a GAP session below, the numbers of non- 
rewritable words exactly coincide with the data from [T]: 

gap> G := AlternatingGroup(5) ; ; 

gap> Exec ("date") ;RewritabilityLength(G, 10) ; time; Exec ("date") ; 

Wed Mar 31 11:04:39 BST 2010 

Started enumeration of NRW of length 2 

29 NRW of length 2 constructed 

Started enumeration of NRW of length 3 

1315 NRW of length 3 constructed 

Started enumeration of NRW of length 4 

43121 NRW of length 4 constructed 

Started enumeration of NRW of length 5 

528069 NRW of length 5 constructed 

Started enumeration of NRW of length 6 

187719 NRW of length 6 constructed 

Started enumeration of NRW of length 7 

1320 NRW of length 7 constructed 

Started enumeration of NRW of length 8 

NRW of length 8 constructed 

8 

33383583 

Wed Mar 31 20:21:09 BST 2010 

Furthermore, it is easy to see that we can process independently each non- 
rewritable word of the length k to derive all non-rewritable words of length bigger 
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than k. Clearly, this allows parallelisation. The parallel version of the algorithm 
was implemented using the master-worker skeleton from the GAP package SCSCP 
[3]. To ensure that no data were lost, the output was modified to return the to- 
tal number of non-rewritable words of each length summing the numbers over all 
parallel procedure calls: 

gap> RewritabilityParallel(AlternatingGroup(5) ,10,4) ; 
[ 0, 29, 1315, 43121, 528069, 187719, 1320, ] 

The third parameter specifies that parallel computation will be started from the 
words of length four. Thus, in the beginning 1315 non-rewritable words of length 
are computed sequentially to ensure an optimal task granularity. 

The program was tested first on an 8-core Intel server, with dual quad-core 
Intel Xeon 5570 2.93GHz / RAM 48 GB / CentOS Linux 5.3 with 1 master and 8 
workers, and then on a cluster consisting of three machines of the same configuration 
as above with 1 master and 24 workers. The average runtime on six measurements 
is given in the table below. 



Number 


Runtime 


Speedup 


Efficiency 


of workers 






= speedup/nr. workers 


4 


3h Om 5s 


3.23 


0.81 


8 


Ih 50m 23s 


5.26 


0.66 


16 


55m 42s 


10.43 


0.65 


24 


37m 16s 


15.59 


0.65 



To summarise, we have provided an independent verification of the result from 
[I] . This may be considered as an additional motivation to find a theoretical proof 
that A^ is 8-rewritable. While the original publication 20 years mentions the us- 
age of PASCAL or C++ for sequential and parallel computations respectively, now 
this computation has been implemented in GAP, being compatible with other GAP 
code and readable by a suitably qualified GAP user. Parallel tools offered in the 
SCSCP package made it possible to parallelise the code within the GAP system 
without the necessity to switch to other traditional for the high-performance com- 
puting languages that support parallelism. Note that an ongoing HPC-GAP project 
(http://www-circa.mcs.st-and.ac.uk/hpcgap.php) is aimed to reengineer the 
GAP system to provide better support for shared and distributed memory pro- 
gramming models, so in the future this example may be hopefully even better 
reimplemented in a new version of the GAP system. 
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Appendix: GAP Source Code for the sequential version 

RewritabilityLength : =f unction (G , limit) 
local eltsG, s, A, orbsA, x, q, n, isnrw, 

nrw, S, i, j, u, K, y, orbsK, v, tw; 
eltsG:=Filtered( G, s -> s <> () ); 
A:=AutomorphismGroup(G) ; 
orbsA:=Orbits(A,G) ; 

x:=Filtered( List(orbsA, q -> q[l] ), q -> q <> () ); 
x:=List( X, q -> [ [ q ], Stabilizer( A,q ) ] ); 
n:=l; 
repeat 

n:=n+l ; nrw:= [] ; 

Print ("Started enumeration of NRW of length ", n, "\n"); 

S := SymmetricGroupC n ); 

S := FilteredC S, s -> s <> () ); 

for i in [1 . .Length (x)] do 

# Print ( i, "/", Length (nrw) , "\r"); 
u := x[i] [1]; K := x[i] [2]; 
if Size(K) = 1 then 

y := eltsG; 
else 

orbsK:=Drbits(K,G) ; 

y:=Filtered( List(orbsK, q -> q[l] ), q -> q <> () ); 

fi; 

tw := u; 
for V in y do 
tw [n] : = V ; 
isnrw :=t rue; 
for s in S do 

if Product (tw)=Product (Permuted (tw,s) ) then 
isnrw := false; break; 

fi; 
od; 
if isnrw then 

Add( nrw, [ ShallowCopy(tw) , 

Intersection( K, Stabilizer( A,v ) ) ] ); 

fi; 
od; 
od; 

Print( Length (nrw) , " NRW of length ", n, " constructed\n") ; 
if nrw= [] then return n; fi; 
X := ShallowCopy( nrw ); 
until n=limit; 
return fail; 
end; 
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